How to project qubits faster using quantum feedback 
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When one performs a continuous measurement, whether on a classical or quantum system, the 
measurement provides a certain average rate at which one becomes certain about the state of the 
system. For a quantum system this is an average rate at which the system is projected onto a pure 
state. We show that for a standard kind of continuous measurement, for a qubit this rate may be 
increased by applying unitary operations during the measurement (that is, by using Hamiltonian 
feedback), in contrast to the equivalent measurement on a classical bit, where reversible operations 
cannot be used to enhance the rate of entropy reduction. We determine the optimal feedback 
algorithm and discuss the Hamiltonian resources required. 

PACS numbers: 03.67.-a, 03.65.Ta, 02.50.Tt, 02.50.Ey 



It was discovered recently (0, that the average 
amount by which the quantum state of a system is pu- 
rified during a measurement depends, in general, on the 
basis in which the measurement is made. Since changing 
the basis of a measurement is equivalent to performing 
a unitary transformation on the system, we can re-state 
this property by saying that a unitary transformation 
may be applied to a system, so as to increase the average 
amount of information that the measurement provides 
about the final (post-measurement) state. This restate- 
ment makes it particularly clear what is being changed 
about the measurement in order to enhance the infor- 
mation: merely the addition of a separate Hamiltonian 
evolution. Taking this point of view, the ability to per- 
form unitary transformations, or equivalently, Hamilto- 
nian evolution, can be considered as a resource which can 
be used to enhance the properties of a fixed measurement 
process. 

Consider now a continuous measurement of, for exam- 
ple, a qubit. This may be represented by a sequence of 
identical 'finite strength' measurements, each of which 
partially projects the qubit onto the basis {|0),|1)} pj 
(We will refer to this as the computational basis, by which 
is simply meant the basis in which information is en- 
coded). As more measurements are made, the purity 
of the state of the qubit increases, until eventually the 
qubit ends up in one of the basis states. Examining this 
process, we find that the state that results from a mea- 
surement in the sequence is not ideal for the purposes of 
purification for the next measurement in the sequence. 
It is therefore possible to perform a unitary transforma- 
tion at the end of each measurement (and depending on 
the measurement result) to increase the average rate at 
which the state is purified. In the continuum limit this 
becomes a continuous feedback process. 

This fact, while interesting as it provides a technique 
for enhancing a continuous measurement, is, if anything, 
more interesting from a fundamental point of view, due 
to the fact that the same measurement process performed 
on a classical bit cannot be enhanced by a reversible 
transformation, since the necessary sequence of transfor- 
mations requires that superposition states of the compu- 
tational basis must be available. Thus, this constitutes 



an example of something which the quantum nature of 
an object makes possible. 

Since this feedback algorithm is designed to enhance 
the properties of the measurement itself, the resulting 
measurement process is an example of an adaptive mea- 
surement, a term first introduced, we believe, by Wise- 
man Q. Wiseman's adaptive measurement scheme, real- 
ized recently by Armen et. al. |^ , was designed to make a 
canonical phase measurement optimally well, given that, 
in practice, one only has arbitrary quadrature measure- 
ments at ones disposal. The topic discussed here is there- 
fore a different kind of application for adaptive measure- 
ments. 

Before we begin the development of the feedback algo- 
rithm, let us say a few words about quantum and classical 
measurements, and the relationship between them. Clas- 
sical measurements on classical systems are described by 
Bayesian inference, and may be written in the same form 
as quantum measurements. If one writes the classical 
probability distribution of the quantity one is measuring 
as a diagonal density matrix, then classical measurements 
are in fact a subset of quantum measurements: While 
quantum measurements are described by a set of opera- 
tors Q n , where the only restriction is that &n&n = b 
classical measurements have the further restriction that 
all the fl n commute with the density matrix describing 
the classical system (for a fuller discussion see [7]). 

The behavior of a quantum measurement with com- 
muting operators therefore reduces to that of a classi- 
cal measurement when either unitary transformations are 
not available to rotate the qubit out of the computational 
basis, or there is a rapid decoherence process which im- 
mediately decoheres the qubit in the computational basis. 
In that case the unitary transformation becomes merely 
a classical diffusion process acting on the bit. 

Here we will be interested in a widely applicable 
model of continuous measurement, where the measure- 
ment record is a Wiener process. We will measure in 
the z-basis of a spin-half system (that is, measure the 
observable a z ), denoting the basis states as {|±)} (thus 
the final result of the measurement will either be |— ) 
or |+)). (Any further reference to 'the computational 
basis' will refer to these states.) To describe the con- 
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tinuous measurement one applies a POVM in each small 
time interval At, where the POVM is chosen to scale 
with time in such a way that a sensible continuum limit 
exists |8J. The resulting continuous measurement is not 
merely a mathematical curiosity, as it corresponds to real 
measurements on physical systems El, 11(1 ITU H2I IT^ | (in 
particular, Korotkov ([H ED, 13 El) 

gives explicit ex- 
amples for solid-state qubits). While the approach in Q 
uses a POVM in each time step which has an infinite 
number of outcomes, for measurements on a two-state 
system one can alternatively employ a POVM with two- 
outcomes. A two-outcome measurement which provides 
information about the computational basis is 

n ± = ^|±)(±| + VT^|t)(tI 



= - [(v 7 ^ + VT^)I ± (V^ - a/I - ~K)a z ] , (1) 

We will refer to this measurement in what follows as 
M.. When k = 1/2 M provides no information about 
the quantum system, leaving the state of knowledge un- 
changed, and when k = or 1 the operators Cl± are rank 
one projectors, so that M. provides the maximal amount 
of information about the final state; this case is an 'infi- 
nite strength' measurement by the terminology of . 

Setting k = 1/2 — y/2j At, and taking the limit of re- 
peated measurements as At — > 0, the evolution of the 
density matrix describing the observers state of knowl- 
edge, p, is given by the stochastic Schrodinger equation 



dp = -7[o-„ [a z ,p]]dt+^2j({a z ,p}+ - 2(a z ))dW, (2) 

where dW is the Gaussian stochastic Wiener increment 
satisfying dW 2 — dt. 

We are interested here in how the observer's uncer- 
tainty of the quantum state reduces over time. In or- 
der to keep the calculations tractable, we will use as our 
measure of uncertainty the so called "linear entropy" , 
P(p) = 1 — Tr[p 2 ]. (That this is a useful measure of un- 
certainty is due to its concavity |2j.) For a completely 
mixed state, P{p) = 1/2, and for a pure state, P(p) = 0. 

In general, the reduction in the linear Entropy will de- 
pend upon the outcome of the measurement. As such, 
a sensible measure of the purifying power of the mea- 
surement is therefore the average reduction in the linear 
entropy, over the two outcomes. From 0] we know that 
the amount by which, on average, the two-outcome mea- 
surement purifies the state (that is, reduces P{p)) is 



(2k- 1)P 



1 - (1 - 2P) cos 2 



1 - (1 - 2P)(2k- l) 2 cos 2 



(3) 



where we have written the initial density matrix as 
p = (/+a-er)/2, with 8 being the angle between a and the 
z-axis, and |a| 2 = (1 — 2P). The reduction in uncertainty 
is maximized when 9 — tt/2. That is, when the basis in 
which the density matrix is diagonal is maximally differ- 
ent from the z-basis, being the basis in which the mea- 
surement is made, the measurement is most effective in 



purifying the state. This somewhat curious result means 
that, if the observer's state-of knowledge is not maxi- 
mally complementary to the measurement basis, and the 
observer's objective is purification, then a unitary trans- 
formation should be applied to the state prior to the mea- 
surement. It is also worth noting that, when the average 
reduction in the entropy is greatest, this reduction is the 
same for both measurement outcomes. Thus, in this case, 
the entropy reduction is deterministic, not random. 

Now let us consider the consequences of this for a 
continuous measurement, performed on an initially com- 
pletely mixed state. Approximating the continuous mea- 
surement by a sequence of measurements M. , we see that 
when we make the first measurement, Ap is maximal. 
However, the result of the first measurement does not 
produce a state diagonal in a basis complementary to 
the z basis, and this is generally true of the state result- 
ing from a measurement in the z-basis. As a result, the 
average purification will not be optimal for at least the 
majority of the measurements M., and thus for essen- 
tially all of the duration of the continuous measurement 
process. 

This suggests, therefore, the following procedure: Af- 
ter each measurement we apply a unitary transformation 
to rotate the state appropriately, so as to achieve max- 
imal Ap for each measurement in the sequence (each 
measurement step). Taking the continuum limit, this re- 
sults in a continuous feedback algorithm which increases 
the rate of projection during the continuous measure- 
ment. Calculating the behavior of the linear entropy, P, 
as a function of time is straightforward, because, as men- 
tioned above, when Ap is maximal, the change in P is 
the same for both outcomes of the measurement. Hence, 
for a finite sequence of steps we have, for the n th step, 



P„ 



(l-b 2 /n) n P , 



(4) 



where b 2 = (2k— 1). In the continuum limit this becomes 



P(t) 



-87* 



P(0). 



(5) 



Without the feedback process, the evolution of Pit) de- 
pends on the measurement outcomes, and as a result is 
stochastic. Thus, not only does the Hamiltonian feed- 
back algorithm increase the rate at which the state is 
purified, but also changes the entropy from being a ran- 
dom function of time to a deterministic one. (In fact, 
it is clear that this is true for any entropy, linear, von 
Neumann, or otherwise.) 

In the absence of feedback, due to its stochastic nature, 
the behavior of Pit) is much harder to obtain, even for 
the simple canonical continuous (classical!) measurement 
we have here. Using the technique in |l4|. one obtains, 
for an initially completely mixed state, 
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cosh(V8fca; 



-dx, 



(6) 



where we use the subscript c to denote the fact that this 
is also the result for a classical continuous measurement 
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FIG. 1: The speed-up factor provided by the optimal Hamil- 
tonian feedback algorithm, as a function of final purity, for a 
continuous measurement in a given basis of a qubit. The case 
displayed is for an initially completely mixed state. 



on a classical bit. An analytic solution for the integral in 
Eq.lJHJ does not appear to exist, and we therefore evaluate 
it numerically. In figure ^ we plot the speed-up factor 
provided by the feedback algorithm in obtaining a given 
final level of purity, when the initial state is completely 
mixed. This factor is independent of 7, and increases 
with the final purity, tending to 2 as the final entropy 
tends to zero. 

The Hamiltonian feedback algorithm above has been 
obtained simply by optimizing the increase in purity of 
the final state at each measurement step. We now show 
that this is, in fact, the optimal feedback algorithm for 
obtaining the maximum purity of the state at any par- 
ticular future time. First, let us denote the map which 
takes us from an initial linear entropy P, to the final lin- 
ear entropy P' , for a single time step, when we use the 
optimal unitary transformation, as f(P). Thus we have 
P' = f(P) = (1 - b 2 )P, so / is linear in P. 

Now, let us consider a single measurement step M, 
where the initial entropy is P. If we use a sub-optimal 
procedure, then one of two states results, and we can 
label the linear entropy of these as Pi = gi (P) and Pi = 
(72(f), respectively. On the other hand, the entropy that 
results from the optimal procedure is f(P), and the three 
entropies satisfy 



> f(P), 



(7) 



where p\ and P2 are the respective probabilities that the 
two sub-optimal entropies were obtained. If we apply the 
optimal procedure to both results (both sides) then since 
f(P) is linear, we have 



£>/(<7i(P)) > fif(P))- 



(8) 



Alternatively, if we were to apply a non-optimal proce- 
dure to both cases (both sides of Eq.Q), then we would 
have, for each of the outcomes i = 1,2, 



J2P^(9i(P)) > f(9i(P))- 
3=1 



(9) 



Thus, by the end of the second step, two sequential non- 
optimal procedures would give 



2 

Yl Pi 

i=l 



Y,PW(9i(P) 



> Y,Pif(9i(P)) 

i=l 

> f(f(P)) 



(10) 



where the first line uses Eq.@, and the second uses the 
linearity of f(P). Thus, after two steps the result of 
non-optimal measurements in both steps gives an average 
entropy which is higher than the result of two optimal 
steps. Clearly this procedure can be repeated n times to 
obtain the result for n measurement steps. Thus, we can 
write the final result of n non-optimal steps as 



> f o{n) (P) 



(11) 



for some pi and Pi. In order for the procedure which 
uses all optimal steps to render a greater entropy than 
a procedure that uses some non-optimal steps, then it 
would have to be possible to apply an optimal step to 
the left hand side of Eq. lfTT)) . such that Y^ =1 Pif{Pi) < 
/°(™+ 1 )(P). However, since f(P) is linear, this is impos- 
sible. We can therefore conclude that the Hamiltonian 
feedback algorithm presented above gives the maximum 
possible entropy reduction for any number of steps, or, 
in the continuum limit, the maximum possible entropy 
reduction for a measurement of any given duration. 

So far we have discussed the effects of the feedback 
algorithm, but not given explicitly the algorithm itself. 
The unitary transformation required after each measure- 
ment step must be such so as to rotate the qubit so that 
the Bloch vector lies in the rcy-plane. The first thing 
to note is that the minimum angle of rotation required 
to do this is achieved by rotating such that the x and 
y elements of the Bloch vector remain in the same pro- 
portions (i.e. by keeping the angle that the Bloch vector 
makes with the x and y axes the same). The second is 
that an application of the measurement M. also keeps 
these angles the same. Assuming that the initial state 
is either completely mixed, or has been rotated prior to 
the measurement so that the Bloch vector lies in the xy- 
plane (which is the optimal thing to do), then the initial 
a z = 0, and the state after a measurement M. is 



Pf 



\{I + \J\ - b 2 {a x a x + a y Oy) ± ba z ). (12) 



The unitary transformation required to rotate such a 
state so that once again a z — 0, using the minimum 
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angle of rotation, is U = exp[— i(a/2)n ■ a], with 



tan 



±6/ (VI - 2PyJ\ - b 2 



n = (— cos(</>), sin(^), 0), 



(13) 
(14) 



1 — 2P, and denned <j> by a x = |a| cos(0). Since 4> remains 
the same throughout the sequence of measurements, n 
remains unchanged throughout the feedback process, and 
it is merely a that changes at each feedback step. After 
the n th measurement P n is given by Eq.Q), and hence 



tan 



±&/(v /l - 2Po(l - b 2 /n) n y/l - b 2 ) 



Note that, when the initial state is completely mixed, the 
rotation angle a, after the first measurement is always 
7r/2, regardless of the strength of M. (i.e. regardless of 
the value of b). 

In the continuum limit, the feedback angle a(t) be- 
comes 



a(t) = tan 



- ^ dW , P(t)<l/2. 
y/l - 2P(0)e~^ t 



(15) 



Since the Hamiltonian required to rotate through a is 
proportional to a(t)/dt, the Hamiltonian required for 
the feedback is proportional to the measurement noise 
dW/dt. Thus the feedback required to obtain an opti- 
mal projection rate is Wiseman-Milburn type Markovian 
feedback, with the addition of a time dependent factor. 
However, this kind of feedback is strictly speaking an 
idealization valid in the limit of a large Hamiltonian, 
since dW/dt is infinite. In addition, the feedback Hamil- 
tonian diverges for an initially completely mixed state, 



since a(0) = 7f/2. Thus, a real continuous feedback pro- 
cedure will provide a lower rate of projection, depending 
on the available Hamiltonian resources. We note that the 
divergence of the feedback Hamiltonian for P = 1/2 has 
an analogue in Wiseman's adaptive phase measurement; 
in that scheme, for this state, it is the rate of change of 
the phase estimate which diverges. It is interesting that 
for both schemes this divergence is associated with the 
fact that the measurement must break the symmetry of 
the state. 

While the feedback algorithm projects a qubit onto 
a final pure state with maximal speed, it is fairly clear 
from the construction of the process that the operations 
corresponding to the many possible final outcomes will 
be mutually non-orthogonal. As a result, this adaptive 
measurement will almost certainly not provide full in- 
formation regarding the initial preparation of the qubit. 
This is why we refer to the algorithm as projecting, rather 
than measuring the qubit, for the sense of the term "mea- 
surement" contains a certain ambiguity: it could mean 
either obtaining information about the initial prepara- 
tion, useful in classical and quantum communication, or 
the final state, useful in quantum feedback control and 
quantum state preparation. While in this case the quan- 
tum state is projected maximally fast, this is at the ex- 
pense of losing information about the initial preparation. 
This raises the question of whether there is a trade-off 
between speed of projection, and loss of initial informa- 
tion, in the kind of measurements considered here. In 
addition, optimal algorithms for projecting higher dimen- 
sional systems, and optimal rates obtainable with fixed 
Hamiltonian resources, are also open questions for fur- 
ther work. 

The author would like to thank Howard Wiseman for 
helpful comments on the manuscript. 
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